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Abstract 

The 14-3-3S are a family of dimeric evolutionary conserved pSer/pThr binding proteins that play a key role in multiple 
biological processes by interacting with a plethora of client proteins. Giardia duodenalis is a flagellated protozoan that 
affects millions of people worldwide causing an acute and chronic diarrheal disease. The single giardial 14-3-3 isoform (gl4- 
3-3), unique in the 14-3-3 family, needs the constitutive phosphorylation of Thr214 and the polyglycylation of its C-terminus 
to be fully functional in vivo. Alteration of the phosphorylation and polyglycylation status affects the parasite differentiation 
into the cyst stage. To further investigate the role of these post-translational modifications, the crystal structure of the gl4- 
3-3 was solved in the unmodified apo form. Oligomers of gl4-3-3 were observed due to domain swapping events at the 
protein C-terminus. The formation of filaments was supported by TEM. Mutational analysis, in combination with native PAGE 
and chemical cross-linking, proved that polyglycylation prevents oligomerization. In silico phosphorylation and molecular 
dynamics simulations supported a structural role for the phosphorylation of Thr214 in promoting target binding. Our 
findings highlight unique structural features of gl 4-3-3 opening novel perspectives on the evolutionary history of this 
protein family and envisaging the possibility to develop anti-giardial drugs targeting g 14-3-3. 
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Introduction 

Eukaryotic 14-3-3s are a family of dimeric highly conserved 
proteins (~30 kDa) with pSer/pThr binding property. When 14- 
3-3s are present in multiple isoforms (e.g. seven in Homo sapiens) 
they can form either homo- and/or heterodimers. Mammals, 
including Homo sapiens, contain seven 14-3-3 genes (i.e., P, y, s, a, 

T, r|), and 15 genes coding for 14-3-3s have been identified in 
the plant Arabidopsis thaliana. [1-2]. 14-3-3s bind to hundreds of 
client proteins thus affecting many cellular processes, including 
metabolism, signal transduction, cell cycle, differentiation, apop- 
tosis, protein trafficking, transcription, stress responses and 
malignant transformation [3]. The overall 14-3-3 architecture is 
well conserved among the different isoforms and eukaryotic 
kingdoms. Each 14-3-3 monomer consists of nine antiparallel a- 
helices (A-I) organized in a cup-like shape with the dimerization 
domain located in the N-terminal half [4—8]. The 14-3-3-target 
interaction is mediated by conserved residues in an amphipathic 
groove of each monomer (formed by a-helices C, E, G, and I) and 
requires specific pSer/ pThr containing motifs on the targets. The 
phosphate moiety of the target binding peptide directiy contacts a 



cluster of positively charged residues in the 14-3-3 amphipathic 
groove including Lys49, Arg56, Argl27 and Tyrl28 (human 14-3- 
3^ numbering), whereas the phosphopeptide backbone and side 
chains establish interactions with a hydrophobic patch of 14-3-3 
conserved residues. Three general 14-3-3 consensus binding motifs 
have been identified: RSX(pS/pT)XP, called mode-1, where pS/ 
pT are phosphorylated serine or threonine residues and X stands 
for any type of residue; RX(Y/F)X(pS/pT)XP (mode-2) [9]; and 
Xp(S/T)Xi_2-COOH (mode-3), in which the phosphorylated 
residue occupies the penultimate C-terminal position in the target 
protein [10]. Nevertheless, 14-3-3s may also interact with some 
non-phosphorylated peptides, such as R18, containing the 
WLDLE sequence [11]. The overall rigid structure of the 14-3-3 
dimer can induce, upon binding, conformational changes in the 
target protein resulting in the inter- and/or intra-compartmental 
sequestration of the target itself, activation/inactivation of the 
target enzymatic activity, and, in few cases, promotion/inhibition 
of the interaction between the target and other proteins [8] . From 
an evolutionary point of view, 14-3-3s can be subdivided into four 
major groups coinciding with the four kingdoms. The existence of 
an ancestral 14-3-3 that evolved independendy into different 
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isoforms in each kingdom has been proposed. In particular, the 
independent clustering of Metazoan epsUon isoforms strongly 
suggesting that the epsUon isoforms are relatively conserved and 
similar to the original animal 14-3-3 ancestor [1]. 

Giardia duodenalis (syn. G. lamblia or G. intestinalis, from now on 
referred to as Giardia) is a flagellated protozoan that parasitizes the 
upper part of the small intestine of mammals, including human, 
and causes giardiasis, the most common parasitic diarrheal disease 
worldwide [12]. Giardia is a deeply branched eukaryote, closer to 
fungi and animals than Euglenozoa [13]. The two life stages of 
Giardia, the binucleated trophozoite, which replicates and colonizes 
the host intestine, and the tetranucleated cyst, the infective and 
resistant stage which can survive in the external environment, can 
be reproduced in laboratory. The phylogenetic position of Giardia 
and its minimalistic genomic and cellular organization make this 
parasite a fascinating model to investigate basic cellular processes 
and different aspects of eukaryotes' evolution [14]. 

Giardia posseses a single 14-3-3 isoform (gl4-3-3) showing high 
sequence identity to the 14-3-3s of the epsilon subgroup. In 
previous works we have demonstrated that gl 4-3-3 is a fully 
functional member of the family with a central role in multiple 
biological pathways of Giardia [15-16]. Interestingly, two post- 
translational modifications (PTMs) make gl4-3-3 unique in the 14- 
3-3 protein family: constitutive phosphorylation and polyglycyla- 
tion. More specifically, gl 4-3-3 is constitutively phosphorylated on 
residue Thr214, which is located in the loop between a-helices H 
and I, and it represents the unique known case of 14-3-3 subjected 
to polyglycylation, which consists in the addition of one or more 
glycine residues to the y-carboxyl group of a glutamate [15]. 

Phosphorylation of different 14-3-3s was previously document- 
ed, but it was observed to only occur in certain physiological 
conditions or cell types and involve residues located either in the 
dimer interface, thus favoring 14-3-3 monomerization, or in the ot- 
helrx G and at the extreme C-terminus, thus modulating the 
binding to the targets [17]. The phosphorylation of Thr214 may 
play a role in promoting the binding of gl4-3-3 to its targets. In 
fact, the gl4-3-3 phospho-mimicking mutant T214E, expressed as 
bacterial recombinant protein, displayed an enhanced binding to 
synthetic phosphopeptides reproducing multiple 14-3-3 binding 
motifs in vitro [18]. Moreover, the expression in Giardia of the non 
phosphorylatable T214A mutant behaved as a dominant negative 
leading to an impaired cyst development. 

As for polyglycylation, this unusual PTM takes place at the 
penultimate gl 4-3-3 C-terminal residue, Glu246, and consists in 
the addition of up to 30 consecutive glycines per monomer. The 
length of the polyglycine chain is stage-dependent and decreases 
down to 10 residues during the cyst formation in parallel with a 
partial re-localization of gl4-3-3 to the nuclei [15,18]. Both the in 
vivo expression of the E246A mutation, which disables the gl 4-3-3 
polyglycylation, and the alteration of the polyglycylation/degly- 
cylation enzyme concentration ratio, affects the intracellular 
localization of the protein and the parasite development into cyst 
[18-19]. 

Furthermore, the structural/functional significance of both 
polyglycylation and phosphorylation has been recendy supported 
by the observation that when the gl 4-3-3 was expressed in 

Drosophila melanogaster the protein resulted devoid of both PTMs 
and was unable to complement fly mutants deleted of either the 
endogenous D14-3-3s or die DLeoII (a 14-3-3(; isoform) [20]. 

Due to the relevance of gl 4-3-3 in many Giardia biological 
processes (i.e. cyst formation) and the peculiar need for constitutive 
PTMs for the protein proper activity in vivo, structural information 
is essential to understand how phosphorylation and polyglycyla- 
tion affect the function of gl 4-3-3. This understanding would 



clearly pave the way for the design of antigiardial drugs targeting 
gl 4-3-3 binding activity. 

In order to better clarify the structural and functional role of 
g 14-3-3 PTMs, we studied the behavior of the protein with and 
without modifications using three complementary approaches. 
First, the crystal structure of a gl4-3-3 in the apo form is presented 
and discussed. The structure shows that the protein, in absence of 
both PTMs, oligomerizes through C-terminal domain swapping 
events. The formation of filaments was also confirmed by 
transmission electron microscopy. Second, the combination of 
various mutants with native PAGE and chemical cross-linking 
indicates that the oligomerization process does not occur when a 
polyglycine stretch is added to the protein C-terminus. Finally, 
molecular dynamic simulations support the hypothesis that the 
phosphorylation of Thr214 promotes the binding to the target(s). 

Materials and Methods 

Vectors' construction 

Escherichia coliJM109 and XLlBlue competent cells were used 
for vector manipulation and recombinant proteins expression, 
respectively. Glutathione S-Transferase (GST)-difopein, GST- 
g 14-3-3 (plasmids pl4-X) and GST-T214E expression vectors 
have been described elsewhere [15,18]. The R200K and the 
T208A mutants were obtained by site-directed mutagenesis using 
the pl4-X vector 15 as template, and the designed primers: 
R200Kfor 5'-AGGGCTTGCGAGCTTGCG AAG AAGGCTT- 
TCGAC-3'; R200Krev 5'-GGCCGCGTCGAA.4GCCTTC- 
TTCGCAAGCTCGCA-3'; T208Afor 5'-GCTTTCGACGCG- 
GCCATCGCAGATCTGGACAAG-3'; T208Arev 5'-CGGTC- 
AGCTTGTCCAGATCTGCGATGGCCGCGTC-3 ' (mutated 
triplets are underlined). The reaction was performed as previously 
detailed [15] and according to manufacturer's instruction. The 
obtained plasmids were designed as "pT208A-X" and "pR200K- 
X". To obtain the polyGlO and polyG20 mutants in which the last 
two C-terminal residues were deleted and replaced with a stretch 
of 10 or 20 glycines, respectively, we took advantage of the 
presence of a Kpnl site at position 404 of gl 4-3-3 coding sequence 
and a NotI site in the multicloning site of the pGEX-6Pl vector. A 
KpnI-NotI cassette was PGR amplified from the pl4-X vectorl5 
using the "gl4KpnIfor" primer, 5'-CCCGGr^CCTCGCCGAG- 
TACTCGT-3' (Kpnl restriction site is in italic), in combination 
with die "polyGlOrev", 5'-GGGGCGGCCGC7TJACCTCCAC- 
CTCCTCCACCTCCACCTCCACC CTCGGCATTATCGT- 
CACC-3', or the "polyG20rev" primers, 5'-GGGGCGGCCG- 
CraACCTCCACCTCCACCTCCACCTCCACCTCCACCT- 
CCACCTCCTCCACCTCCACCTCCACC CTCGGCATTAT- 
CGTCACC -3' (sequence coding for polyglycines stretch is in bold 
and underlined, the NotI restriction site is in italic and the stop 
codon is underlined and in italic). The coding sequence of human 
14-3-3i; was PGR amplified from the plasmid pGEX2TK-14-3-3i; 
[21] using the primers "14^forF", 5 ' -C GC G CATC C ATGG A- 
TAAAAATGAGCTGGTTC-3' (BamHI restriction site is under- 
line), and "14Crev", 5'-GGGGCGGCCGCTTAATTTTCCCC- 
TCCTTCTCC-3' (Nod restriction site is in italic). PGR reactions 
were performed in a final volume of 50 |il using 25 |il of 2X PGR 
master mix (Promega, France), 20 pmols of each primer and 50 ng 
of plasmid pl4-X [15] as template. Reactions were performed on 
a T-Personal Thermocycler (Biometra Corporation, Gottingen, 
Germany). Amplification conditions were: one cycle at 95°C for 
5 min; 30 cycles at 95°C for 30 sec, 55°C for 30 sec and 72°G for 
30 sec; and one cycle at 72°G for 7 min. The PGR fragments were 
firstiy cloned in the pGEM_Teasy vector (Promega). KpnI/NotI 
digested polyGlO and polyG20 fragments were sub-cloned in the 
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KpnI/Notl-digested pl4-X vector replacing the last 341 nucleo- 
tide of the gl4-3-3 coding sequence, whereas the BamHI/NotI 
digested 14-3-3^ was cloned in the BamHI/Notl-digested 
pGEX6Pl vector. The obtained plasmids were designated as 
"pPolyGlO-X", "pPolyG20-X" and "phl4-3-3i;-X", respectively. 

Proteins expression and purification 

For protein expression, transformed E. coli were grown in SOB 
medium at ODgoo = 0.6-0.8, and expression induced with 0.,5 mM 
isopropyl diio-|3-D-galactoside (IPTG) at 37°C for 4 h. GST-fused 
proteins were purified by affinity chromatography on glutathione- 
sepharose 4B (GE Healthcare, Uppsala, Sweden) and released 
from GST by digestion with PreScission protease (GE Healthcare) 
at 4°C for 16 h in digestion buffer according to manufacturer. 
Protein were diafiltrated o.n. at 4°C either in 50 mM Tris-HCl 
pH 7.5 (for crystallization) or in 10 mM Tris-HCl pH 7.0, for 
Circular Dicroism (CD) experiments, using a PM-5 membrane. 
Proteins were concentrated using Centricon 10 (MiUipore 
Corporation, Bedford, MA, USA) and concentration was 
measured with Bradford's method (BioRad, Hercules, CA, USA). 

Circular dicroism 

CD spectra were measured with a Jast:o J-715 spectropolarim- 
eter (Jasco Ltd, Hachioji City, Tokyo, Japan) in 1.0 cm quartz 
cuvettes between 260 and 195 nm. Samples were 0.1 Itg/^l in 
10 mM Tris/HCl buffer (pH 7.0). 

Protein crystallization, data collection and processing 

Crystals were grown by the hanging drop vapor diffusion 
method at 293 K. Protein sample was concentrated to 10.5 mg/ 
ml. Crystallization drops consisted of 1 |tl of protein solution 
mixed with an equal volume of the reservoir solution on a cover 
slip which was suspended over the reservoir containing 2 M 
ammonium sulfate, 5"/, (w/v) PEG400 and 0. 1 M HEPES pH 7.5. 
Crystals, grown in 1-2 weeks, were cryo-protected with 20% 
glycerol, mounted in nylon loops and flash-cooled by submersion 
into liquid for transport to the synchrotron-radiation source. 
Crystals are poorly reproducible and most of them diffract at low 
resolution. X-ray diffraction data were coUected at 3.2 A 
resolution as 0.5° oscillation frames at 100 K on the beam line 
ID 14-1 at ESRF (Grenoble, France) using a CCD detector. The 
data were processed using the HKL2000 package [22] and the 
CCP4 suite [23]. The structure was solved by molecular 
replacement with the program Phaser [24] using the structure of 
human 14-3-3E as search model (PDB code 2BR9) [25] which 
displays 63% sequence identity with gl4-3-3. Manual fitting and 
model building were performed using COOT [26]. Refinement 
was performed using Refmac5 [27]. Structural figures were 
generated with PyMol [28]. Crystsd parameters, data-coUection 
and refinement statistics are presented in Table 1. Atomic 
coordinates and structure factors have been deposited in the 
Protein Data Bank with accession number 4F7R. Protein 
Interfaces, Surfaces and Assemblies (PISA) Service [29] at the 
European Bioinformatics Institute (http://www.ebi.ac.uk/msd- 
srv/prot_int/pistart.html) was used for protein-protein complex 
analyses. The structural alignment between proteins was per- 
formed using the software Superpose (Secondary Structure 
Matching) [30]. 

Molecular dynamic simulations 

The in silico phosphorylation of Thr2 1 4 was modeled on the 
three-dimensional crystal structure of gl4-3-3 (PDB:4F7R) using 
the psfgen module of NAMD [31]. The molecular dynamics 



simulation was performed using the NAMD 2.8 package for GPU 
computing [31]. The gl 4-3-3 structure was immersed in a 
rectangular simulative box fiHed with 45.000 TIP3P [32] explicit 
solvent molecules and rendered electro neutral by the introduction 
of 4 sodium counter-ions. A standard CHARMM force field [33] 
was adopted to parameterize the system for bond parameters and 
Lennard-Jones terms [34]. The length of hydrogen bonds was 
constraint by applying the LINCS algorithm [35] and the 
electrostatic interactions were taken into account by means of 
the Particle Mesh Ewald method [36] . A first round of 5000-step 
minimization was performed in order to regularize the structure of 
aU the atoms composing the system, using the Steepest Descent 
algorithm [37]. Optimization and relaxation of solvent and ions 
was performed by simulating the system at a temperature of 50, 
100, 150, 200, 250 K, for 500 ps each with a time step of 1.0 fs. 
Simulations were performed at a constant temperature of 300 K 
using the Langevin method and at constant pressure; the pressure 
was kept constant (1 bar) using the Rahman-ParrineUo barostat 
[38] and 2.0 fs time step. Principal component analysis (PCA), 
structural clustering, and analyses of RMSD, RMSF, solvent 
accessible surface area, and salt bridges were carried out using the 
GROMACS 4.5.4 analysis tools" [39]. PCA (i.e. Essential 
Dynamics) was performed as previously described [40]. Briefly, 
higher frequency fluctuations were filtered out by diagonalizing 
the atomic positional fluctuation covariance matrix calculated 
during the production run. The corresponding eigenvectors and 
eigenvalues were used to describe large-amplitude motions. The 
first few eigenvectors of the diagonalized covariance matrix usually 
account for the major fraction of the total variance. The projection 
of the atomic trajectories onto the corresponding eigenvectors 
represents large collective atomic motions. The images were 
obtained with PyMol and UCSF Chimera [41] and the graphs 
were calculated with XMGRACE (Paul J. Turner Center for 
Coastal and Land-Margin Research Oregon Graduate Institute of 
Science and Technology Beaverton, Oregon). 

Purification of endogenous gl4-3-3 

Trophozoite of Giardia duodenalis WBC6 strain were grown as 
previously described [15]. Parasite's protein lysate was prepared as 
previously described [18] and endogenous g 14-3-3 was affinity 
purified on glutathione-Sepharose immobilized GST-difopein and 
gl4-3-3 eluited with 100 |tl of 5 mM "Raflp" phosphopeptide, 
(LSQRQRST(pS)TPNVHMV), reproducing one of the 14-3-3- 
binding motifs of the human Raf-1 protein [15]. The excess of 
phosphopeptide was removed by dialysis against HT buffer 
(25 mM Hepes-KOH pH 7.6, 75 mM KCl, 5 mM MgClg, 
0.1 mM EDTA, 0.04% Tween 20) with Centricon 10 (Millipore 
Corporation, Bedford, MA, USA). 

Mass spectrometry analysis 

Affinity purified gl4-3-3 was separated on ID-gel NuPAGE 4— 
1 2 % (Novex, Invitrogen) run in morpholinepropanesulfoninic acid 
(MOPS) buffer and stained with the Colloidal Blue Staining kit 
(Invitrogen). For the MS/MS identification of phosphorylated 
residue in gl433, in gel digestion and nanoflow reversed-phase 
liquid chromatography tandem mass spectrometry (RP-LC-MS/ 
MS) analysis were performed, with minimal changes, as previously 
described [16] using an HPLC Ultimate 3000 (DIONEX, 
Sunnyvale, CA U.S.A) on line connected with a linear Ion Trap 
(LTQ,-XL, ThermoElectron, San Jose, CA). MS data were 
searched for matching with the Giardia protein database, 
downloaded from the web site http://www.giardiadb.org/ 
giardiadb, considering peptides specificaUy cleaved by trypsin on 
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Table 1. Crystal parameters, data collection statistics and refinement statistics of g14-3-3. 




PDB ID 


4F7R 


Space group 


P32 


Unit cell parameters (A) 




1 QQ g 




1 4Q 5 


No. of molecules in the asymmetric unit 


4 


IfJt dLUIIIK. MIULICl \r\ j 


54 5 


Resolution ranges (A) 


'^n n ^ 7 n ^1 ^ 9ni 

JU.U J.Z ^J.J 1 D.£.\J) 


Unicjue reflections 


26426 


Completeness (%) 


1 no ^^ nr\\ 


Redundancy 


7.0 (7.0) 


"merge 


12.5 (76.2) 




1.06 (1.2) 


<//a(/)> 


16.3 (2.8) 


Rcrys (%) 


18.5 (29.4) 


Rfree (%) 


23.1 (31.9) 


rms (angles) (°) 


0.812 


rms {bonds} (A) 


0.005 


Residues in core region of Ramachandran plot (%) 


930 (98.9) 


Residues in generously allowed region of Ramachandran plot (%) 


10 (1.1) 


Residues in disallowed region of Ramachandran plot (%) 


0 


Values in parentheses are for the highest-resolution shell. 

""merge = Z),hE/I'/(''W) - </(/iW)>|/E/.«E'«'''<"' '"'^^'^ is the /th observation of the reflection [hkl] and <l(hkl)> 
Y = h,{\lij(hkD - <mD>\f/{ai'N/{N -1). 
doi:1 0.1 371 /journal.pone.0092902.t001 


is the mean intensity of the [hkl) reflection. 



K or R, possibly modified by oxidation on methionine and/or 
phosphorylation on S, T or Y residues (Am: -I- 80 Da). 

Cross-linking Studies and native page 

For cross-linking studies, recombinant proteins were dialyzed in 
HT bufTer with Centricon 10 (MUlipore Corporation, Bedford, 
MA, USA). A final concentration of 0.33 |iM recombinant 
proteins, or affinity purified endogenous gl4-3-3, were pre- 
incubate with "Raflp" phosphopeptide (0.66 |iM) in 5 |il of HT 
bulfer for 16 h at 4°C either in presence or absence of 1 mM DTT 
and then subjected to cross-linking by the addition of 5 |il of 0,4 M 
triethanolammine pH 8.0 with or without 300 DMP 
(dimethyl pimelimidate) for 5 min at 37°C. Reaction was stopped 
with 5 |il of 150 mM Tris-HCl pH 8.0. Samples were mixed with 
5 |ll of 4X loading buffer (Novex, Invitrogen) and 2 |il of sample 
reducing agent (Novex, Invitrogen), separated on NuPAGE 4— 
12% gel (Novex, Invitrogen) in morpholinepropanesulfonic acid 
(MOPS) buffer, transferred on nitrocellulose (BioRad, Hercules, 
CA, USA) and probed with the rabbit anti-gl 4-3-3 serum (serum 
N14 at 1:5000) [15] or the rabbit anti-human 14-3-3 antibody 
Abl4112 (at 1:500, Abeam pic, Cambridge, UK). Antibody/ 
antigen interaction was revealed by incubation with anti-rabbit 
HRP-conjugated secondary antibody (Ab) (1:2,000) and chemio- 
luminescence (MiUipore Corporation, Bedford, MA, USA). For 
native PAGE, 5 |iM of recombinant proteins or affinity purified 
endogenous g 14-3-3 were incubated in HT buffer with 10 [xM of 
"Raflp" peptide or the "Rafl" unphosphorylated one 16 h at 
4°C. Native PAGE was performed on 12% PAGE as previously 
described [42] and gel stained with 0.25% (w/v) Coomassie 
Brilliant Blue R-250. 



Transmission Electron Microscopy (TEM) 

Recombinant Prescission-cleaved gl 4-3-3 and polyG20 mutant 
were incubated at 4°C at 0. 1 or 1 |ig/ |il in 25 mM Hepes-KOH 
pH 7.6 for 16 h. Aliquots of 10 |il of protein solutions (0.1 \ig/[il 
and 1 |.tg/|il) were apphed onto formvar-carbon coated copper 
grids. After 5 minutes for adsorption, the excess of protein solution 
was sucked off The air-dried grids were stained with an aqueous 
solution of 1 % (w/v) uranyl acetate for 15 s. The excess of staining 
solution was removed carefully. The grids were observed with a 
PHILIPS EM208S transmission electron microscope at an 
acceleration voltage of 80 kV. Electron micrographs were taken 
with a slow scan camera (Megaview III, Olympus). 

Accession number 

The PDB accession number for the gl4-3-3 reported in this 
paper is 4F7R. 

Results 

The overall structure of the gl4-3-3 

To achieve an amount of gl 4-3-3 suitable for structural analyses 
we used an N-terminally GST-fused gl 4-3-3 expressed in bacteria 
that was characterized in previous works for its ability to bind both 
to several phospliopeptides reproducing 14-3-3 binding motifs and 
to the phosphorylation-independent binding peptide difopein 
[15,18]. Although the gl4-3-3-specific polyglycylase gTTLL3 
(Tubulin Tyrosin Ligase-Like3) had already been identified and 
expressed [19], the low efficiency of the in vitro polyglycylation 
process did not allow us to to obtain sufficient polyglycylated 
protein for crystallization studies. On the other hand, the gl4-3-3- 
specific protein kinase is still unknown. Notice that the GST- 
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Figure 1. Classical dimer of g14-3-3 assumes an "open" conformation. A) Two views of g14-3-3 N-terminal dimer (green, cyan). B) The 
superimposition of g14-3-3 and human 14-3-3E (2BR9, violet) emphasizes the open conformation of g14-3-3. C) Close up of N-terminal dimerization 
interface. Residues are shown as lines if involved in hydrogen bonds, as sticks if involved in salt bridges. D) Detail of polar interactions between 
helices H and I. Residues not conserved in 14-3-3 family are represented as sticks. A conserved residue (Gln226 in g14-3-3 Gln222 in human 14-3-38) is 
shown as yellow sticks in B) and D) to highlight the different orientations of helix I. 
doi:1 0.1 371 /journal.pone.0092902.g001 



cleaved full-length unmodified gl4-3-3 used for the crystaUization 
experiments only diHered from the native gl 4-3-3 for the presence 
of a GPLGS sequence upstream the initial methionine. 



The structure of the gl4-3-3 was determined at 3.2 A resolution 
in the P32 space group and the asymmetric unit contained two 
functional dimers (Fig. lA-lB). Each gl4-3-3 monomer displayed 
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the typical 14-3-3 structural organization [5-7] with nine 
antiparallel a-helices (A-I) arranged in a 'U' shaped fashion with 
a large central binding groove. Each monomer could be divided 
into two sub-domains: the N-terminal sub-domain (helices A-D), 
which consisted of the canonical dimerization interface as well as 
the floor of the binding groove, and the C-terminal sub-domain 
(hehces E-H), which comprised the side walls and the roof of the 
binding groove (Fig. lA-lB). The gl4-3-3 buried N-terminal 
dimerization interface was 1240 A^ per monomer (8.5% of the 
total accessible surface area), scored 0.3 in Complexation 
Significance Score (CSS), indicating that the dimc-rization was 
not due to crystal packing and confirming, as expected, the 
relevance of the gl4-3-3 dimerization surface. The gl4-3-3 N- 
terminal dimerization relies on residues largely conserved among 
the 14-3-3 family members [6,25] involving, in this case, three salt 
bridges between Arg22 and Glu97, Asp9 and Lys82, and Glu25 
and Lys93 (Fig. lA and IC). Arg22, Asp9, and Glu25 are located 
in gl4-3-3 a-helices A and B of one monomer whereas Glu97, 
Lys82, Lys93 reside in a-hehx D of the other protomer. In the 
homodimeric human 14-3-3S, as well as in the epsUon-like Cpl4e 
from the protozoan parasite Cryptosporidium parvuni, dimerization 
involves only the first salt bridge (Argl9-Glu92 in hl4-3-3s and 
Arg43-Glull9 in epsilon-like Cpl4e, respectively) [7,25]. The 
presence of a single salt bridge allows hl4-3-3s to mainly 
heterodimerize with other human isoforms [25]. Interestingly, 
gl4-3-3 dimerized through three salt bridges similar to 14-3-3^ [4- 
41] despite the higher sequence identity of g 14-3-3 with proteins of 
the 14-3-3S subgroup than with those of the 14-3-3^ subgroup. 
This suggests that gl 4-3-3 only forms stable homodimers. In 
support of this, when g 14-3-3 was expressed in Drosophyla 
melanogaster or when fly isoforms were expressed in Giardia 
heterodimerization occurred in both cases only between gl 4-3-3 
and D14-3-3i; (LeoII) but never between gl4-3-3 and D14-3-3E 
[20]. The three-dimensional superimposition of the gl 4-3-3 and 
h 14-3-38 structures yielded a root mean square deviation (rmsd) of 
0.95 A. This low rmsd value indicates that helices A-H are largely 
super-imposable. The software automatically excluded helix I 
from the calculation as it was recognized as a largely dissimUar 
region. Indeed, helix I of gl4-3-3, compared with the human 
isoforms, is rotated of about 180° around helix H (Fig. IB and ID) 
implying a substantial rearrangement of the loop connecting 
helices H and I (HI loop). The displacement of the HI loop causes 
the disruption of a "ST-motif' found in most 14-3-3s and 
including Ser210 (according to hl4-3-3i^ numbering), which 
corresponds to Thr214 in gl4-3-3. (Fig. ID). This motif consists 
of 4—5 residues including a Ser or a Thr forming hydrogen bonds 
with two residues downstream in the sequence [43]. Instead, the 
displacement of the HI loop in gl4-3-3 is favored by eight 
hydrogen bonds between helices H and I (Fig. ID) that involve, 
among otiiers, residues Arg200, Thr208, Thr214 and Asn233 
which correspond, in the human isoforms, to Lysl94, Ala202, 
Ser/Asn208 and Thr227, respectively (according to hl4-3-3i^ 
numbering) (Fig. 2). 

As consequence of the helix I rearrangement, the gl4-3-3 
structure showed a fuUy open "extended" conformation of the C- 
terminal region, supporting the notion that the g 14-3-3 region 
encompassing helices G and I possesses a degree of flexibflity 
allowing the binding groove to change from a "closed" to an 
"open" conformation and thus accommodating different types of 
peptides [25]. 

C-terminal dimer-dimer interface 

The described gl4-3-3 "extended" conformation entafls the 
exposure of several hydrophobic residues located in helix G, H 



and I and usuafly buried in other 14-3-3s. A possible consequence 
of such exposure of hydrophobic residues is the observed domain 
swapping occurring between the C-terminus of two gl4-3-3 
monomers, each belonging to a different dimer, in which helix I' 
from a monomer interacts with helices F and H from the 
monomer of the neighboring dimer (Fig. 3A; B and C). In contrast, 
in other available structures of 14-3-3 dimers, hehx I is observed to 
interact with helices F and H from the same monomer [6]. 
Domain swapping entails a dimerization event at the protein C- 
terminus with the formation of a four-helix buncUe between two 
adjacent gl 4-3-3 dimers (Fig. 3B and 3C), leading to the 
assemblage of protein filaments. Domain swapping was observed 
in the structure of different proteins and it is defined as two 
identical protein chains exchanging a part of their structure to 
form an intertwined dimer or a higher-order oligomer [44]. The 
overall g 14-3-3 crystal packing resulted in stacks of endless paraUel 
filaments (Fig. 3D), each layer being rotated 60° with respect to the 
previous one. The calculated C-terminal dimer-dimer buried 
interface is 1576.9 per monomer, i.e. about 11% of the total 
accessible surface area. The calctdated solvation free energy 
gained upon formation of the interface (AG — -39.3 kcal/mol) 
suggests that the oligomer is highly stable. As previously stated, a 
stretch of hydrophobic residues is exposed at the C-terminus of 
each monomer so that the C-terminal dimer-dimer interface 
results predominantiy hydrophobic and further stabilized by four 
intermolecular hydrogen bonds between Ser221 of helix I and 
Thrl70 of loop FG' (and Ser221 of helix Y and Thrl70 of loop 
FG), and Glu215 of loop HI and Leu213 of helix H' (and Glu215 
of loop HI' and Leu2 1 3 of helix H). Nearly, all the described 
intermolecular interactions in g 14-3-3 can be also observed in the 
available 14-3-3 structures in the form of intra-molecular 
interactions involving helix I, which perfectiy replaces the swapped 
helix r (Fig. 3B). The only exception is represented by the 
hydrogen bond involving loop HI (residues Glu215-Leu213) that 
does not occur in the non-swapped dimer (Fig. 3B). 

The gl 4-3-3 target binding site 

The 14-3-3 binding to its target(s) primarily occurs through the 
interaction of conserved residues in the 14-3-3 binding groove and 
a phosphorylated peptide in the target(s). According to hl4-3-3E 
numbering, the pSer/pThr phosphate moiety of the ligand binds 
to a highly conserved 14-3-3 triad formed by Arg57, Argl30 and 
TyrlSl whilst Asnl76 and Asn227 interact with the phosphopep- 
tide backbone [6,2.'i]. 

In the g 14-3-3 structure, despite helix I swapping, the 
conformation of the residues involved in the recognition and 
binding to phosphopeptide(s) resulted to be conserved (Fig. 4). 
Indeed, the visual inspection of the Fo-Fc electron density map 
revealed a peak that was assigned to a sulphate ion (Sup. Fig. SI), 
present in the crystallization solution, which lies in the phosphate 
binding site and interacts with the gl4-3-3 Arg60-Argl35-Tyrl36 
triad (Fig. 4A and Sup. Fig. SI), thus mimicking the phosphate 
moiety of a pSer/ pThr binding peptide. This was supported by the 
superimposition between g 14-3-3 (in the apo form) and the human 
14-3-3e/RRQR(pS)AP phosphopeptide complex (PDB 2BR9). All 
the residues involved in the phosphopeptide binding (Arg60, 
Argl35, Aspl32, Tyrl36, Asnl80, Asn231, according to gl4-3-3 
numbering) are conserved, although it is worth noting that Asn23 1 
belongs to helix I' in gl4-3-3 (Fig. 4B-4C). The observed 
"extended" conformation of gl4-3-3 with the large conforma- 
tional change of helix I would disrupt the binding groove, which is 
indeed restored by domain swapping of helix I', indicating that the 
protein assemblies can retain phosphopeptide binding properties 
even in the case of C-terminal dimerization (Fig. 4B— D). 
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Figure 2. Sequence alignment of g14-3-3 protein. The g14-3-3 protein has been aligned with two Ordeum vulgaris isoforms, one from epsilon 
and one from non-epsilon plant 14-3-3 subgroups (b14-3-3A, P29305.1; b14-3-3B, Q43470.1, respectively), the two Saccharomyces cerevisiae yeast 
isoforms (BMHl, P29311.4; BIV1H2, P34730.3) and the seven human 14-3-3 isoforms (h14-3-3|3, accession number NP_647539.1; h14-3-3T, 
NP_006817.1; h14-3-3ri, NP_003396.1; h14-3-3C, NP_003397.1; h14-3-3y, NP_036611.2; h14-3-3a, NP_006133.1; h14-3-3c, NP_006752.1). Alignment 
has been generated using ClustalW2 software and edited with BOXSHADE 3.21. Identical residues are black boxed, similar residues are gray boxed, 
divergent ones are left unboxed. Dashes indicate gaps. A grey line above the alignment indicates the a- helices, the thin line indicates a region of 
helix3,o, observed also in h14-3-3 y, 8, t| and a. Residues involved in N-terminal dimerization are indicated by stars; the triad of residues contacting 
the phosphate moiety in the target phosphopeptide are indicated by empty triangles. Dots indicate other residues taking contact with the target 
phosphopeptide backbone and lateral chains. Black triangles indicate Thr214 and Glu246 of g14-3-3. White squares indicate the residues Arg200, 
Thr208, and Asn233 forming intermolecular hydrogen bound in the C-terminal dimer interface. 
doi:1 0.1 371/journal.pone.0092902.g002 



Structural role of g 14-3-3 phosphorylation 

The solved stmctur(; of gl4-3-3 is not phosphor^'lated and 
Thr214 is buried at the C-terminal interface (Fig. ID) and its side 
chain establishes a hydrogen bond with Thr222. Therefore, it 
seems Ukely that phosphorylation alters the Thr214 pattern of 
interactions (i.e. potentially affecting the local protein folding or 
energy) thus influencing target binding and/ or domain swapping. 
In order to better understand how phosphorylation of Thr214 
may affix! the protein motion, the effects of phosphorylation on 
the dynamical behavior of g 14-3-3 were analyzed by in silico 
phosphorylation of Thr214 in the 3D structure of the gl4-3-3 
monomer and MD simulations on both the unmodified (WT-gl4- 
3-3) and the phosphorylated systems (Ph()-gl4-3-3). An MD 
simulation of 50 ns was run and the first 20 ns were discarded to 
take into account the system equilibration so that all the analyses 
were carried out on the last 30 ns (Fig. 5A). The effects of the 
phosphate group on the protein motion are clearly visible (Fig. 5B). 
In the WT-gl 4-3-3 monomer, helices H and I fluctuated much 
more than in Pho-gl 4-3-3. The calculated per residue Root Mean 
Square Fluctuation (RMSF) showed that residues 1-170 fluctuated 
similarly in WT-gl4-3-3 and Pho-gl4-3-3, whereas differences 



could be observed in the helices upstream and downstream 
Thr214 (Fig. 5B). Interestingly, fluctuations of the loop HI linking 
helices H and I were similar in the two systems. The behavior of 
the C-terminal of the protein upon phosphorylation was even 
more evident in the principal component analysis (PCA) results 
(Fig. 5C, 5D and 5E). The first 10 eigenvectors (Fig. 5B) described 
about the 70% of the fluctuations in the essential space and 
represented most of the protein motion in both phosphorylated 
and non-phosphorylated proteins (each dot in the plot represents 
an eigenvector). However, the weight of the first two eigenvector 
represents about the 35% of the global fluctuations in the WT- 
gl4-3-3 monomer (Fig. 5C), and only the 25% in Pho-gl4-3-3 
(Fig. 5C). The first two eigenvectors described predominantly the 
motion of helices H and I, which were subjected to larger 
movements in WT-gl4-3-3 (Fig. 5E) than in Pho-gl4-3-3 (Fig.SD), 
where they were also closer to each otiuT. The r(;sults of the PCA 
were consistent with the MD snapshots taken at 20, 30, 40, and 
50 ns (Fig. 6A). Overall, the computational study showed that the 
phosphorylated residue acts as a "lock" in the elbow between helix 
H and I. The negatively charged phosphate, indeed, could stiU 
allow the interaction with the polar side chain of Thr222, as the 
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Figure 3. C-terminal domain swapping leads to fibrils forma- 
tion in the crystal. A) Two adjacent dimers within the crystal (green, 
orange), corresponding to one asymmetric unit, give rise to non- 
canonical C-terminal dimerization through the swapping of helix I. B) 
Detail of the interface showing the residues involved in hydrogen 
bonds. C) The superimposition of g14-3-3 and human 14-3-3e (violet) 
emphasizes the swapping of helix I (g14-3-3 otl' vs h14-3-3E otl) as 
indicated by the arrow. D) Crystal packing (36 unit cells, 6x6x1). The 
crystal is formed by layers (blue, red, yellow) of endless parallel 
filaments. 

doi:1 0.1 371 /journal.pone.0092902.g003 

non-phosphorylated Thr214 did, and additionally formed a 
network of interactions with the polar side chain of Tyr2 1 8 and 
with the positively charged side chain of Lys219 (Fig. 6B), by 
estabhshing several hydrogen bonds and a salt bridge. The 
stiffening of the hinge between helices H and I in presence of the 
phosphate group forced the two helices to keep a much closed 
conformation during the simulation and oscillated much less than 
in the non-phosphorylated protein. These interactions were much 
weaker, or even missed, in WT-gl4-3-3 where the negatively 
charged phosphate group on Thr2 1 4 is absent. The lack of such 
interaction network in WT-g 14-3-3 made helix I more flexible and 
capable of moving away from helix H, thus permitting the 
observed "extended" conformation of the protein C-terminal. 
Despite the more closed conformation imposed by the phosphor- 
ylation of Thr214, the structural .superimposition of Pho-g 14-3-3 




Figure 4. Phosphopeptide binding site is preserved upon helix 
I swapping. A) A sulfate ion bound in correspondence of the 
phosphate binding site in g14-3-3. 2mFo-dFc electron density 
(contoured at la) is shown for sulfate ion and the residues involved 
in the interaction. B) Superimposition of g14-3-3 (green, orange) and 
14-3-3s-phosphopeptide (2BR9, violet-pink) complex showing that 
the swapped helix 1' can participate to peptide binding. C) Detail of 
the peptide-binding site. The peptide (from the 2BR9 structure) and the 
residues involved in the interaction, all conserved in g14-3-3, are 
represented as sticks. D) Orientation of the phosphopeptides modeled 
into a g14-3-3 filament. The phosphopeptides in the binding cleft are 
represented by arrows blue to red(N-terminus to C-terminus). 
doi:10.1371/journal.pone.0092902.g004 



(at 20 ns of MD), WT-gl4-3-3 and the hl4-3-3E/RRQRpSAP 
complex (Fig. 6C) shows that the Pho-gl4-3-3 helix I does not 
overlap with helix I of hl4-3-3s, as indeed the swapped helix V of 
g 14-3-3 does (Fig. 4), leading to the possibility that phosphoryla- 
tion alone might not be sufficient to both hide the patch of 
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Figure 5. Molecular Dynamic simulation. A) The evolution of the RMSD (Root Mean Square Deviation) of WT-g14-3-3 (in red) and Pho-g14-3-3 (In 
black) during the simulation. In order to allow the system to equilibrate, the first 20 ns of simulation were discarded from the analyses (region shaded 
in grey). B) The per residue RMSF (Root Mean Square Fluctuation). Helices H and I are enclosed in dashed boxes. The position of Thr214 is indicated 
with an arrow. The bottom plot (grey) shows the per residue RMSF difference between WT-g14-3-3 and Pho-g14-3-3. A schematic secondary structure 
diagram of the gl4-3-3 is reported on the top. C) Eigenvectors cumulative weight on total motion for WT-g14-3-3 (in red) and Pho-g14-3-3 (in black). 
Projection of motion along the first eigenvector for Pho-g14-3-3 D) and the WT-g14-3-3 E) protein structures. The amplitude of motion follows the 
color scale from red to blue. The phosphorylated residue is represented in grey and indicated with an arrow. 
doi:1 0.1 371 /journal.pone.0092902.g005 



hydrophobic residues that are in the C-terminal dimer interface 
and impair the C-terminal domain swapping. 

Effect of amino acid substitutions and PTMs on 
oligomerization 

In order to confirm tlie tendency of recombinant gl 4-3-3 to 
form multimers and shed light on the mechanisms leading to the 
observed "extended" conformation and domain swapping, differ- 
ent gl 4-3-3 mutants were constructed. Mutants were analyzed by 
native PAGE and chemical cross-linking, two well established 
methods already used to asses hi 4-3-3^ dimerization properties 
[42]. As previously observed, residues Arg200, Thr208, and 
Asn233, which are unique to gl4-3-3, may contribute to promote/ 
stabilize the "extended" conformation and the co-occurring C- 
terminal domain swapping. For this reason, two of them were 
mutated, Arg200 into lysine (R200K) and Thr208 into alanine 
(T208A). Lysine and alanine are the residues found in the 
corresponding positions in all human 14-3-3 isoforms (Fig. 2). 
Furthermore, due to the relevance of polyglycylation, two mutants 
were constructed in which the last two C-terminal residues of g 1 4- 
3-3 were replaced by a stretch of 10 or 20 glycines (polyGlO and 



polyG20). These mutants could potentially mimic both the 
presence and the variation in length of the polyglycine chain 
occurring during Giardia encystation [15]. In addition, the 
previously described T214E mutant was also included [18]. As 
controls, the well-characterized hi 4-3-3^ and the affmity-purified 
endogenous gl4-3-3 from Giardia trophozoites, which are physi- 
ologically phosphorylated at Thr214 (Sup. Fig. S2), and with 
different length of polyglycine chain at Glu246 [15], were also 
analyzed. No variation in the secondary structure could be 
evidenced by circular dicroism analyses between the unmodified 
recombinant gl4-3-3 and the mutant proteins (Fig. 7A). 

As a result of the chemical cross-linking reactions using DMP 
(9.2 A spacer arm), bands compatible with the formation of 
covalent dimer products were visible for both the endogenous gl4- 
3-3 and, as expected, for the hl4-3-3i^ [42], as consequence of the 
N-terminal dimerization (Fig. 7B, left panel). In agreement with 
the solved crystal structure, dimeric, trimeric and tetrameric forms 
were clearly visible in the recombinant g 14-3-3, and each of the 
T214E and R200K mutants, likely as combination of both N- 
terminal and C-terminal dimerization. Moreover, faint bands at 
molecular weight compatible with pentameric forms were also 
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Figure 6. Molecular Dynamic simulation. A) Snapshots at 20, 30, 40, and 50 ns of the protein structures during the simulation. The proteins are 
superimposed to each other, WT-gl4-3-3 is in green and Pho-gl4-3-3 is in violet. B) Magnification of the loop between helices H and I. The 
phosphorylated Thr214 and the surrounding polar residues are represented in sticks colored by atom. C) Superimposition Pho-g 14-3-3 at 20 ns of IVID 
(violet), WT-g14-3-3 (green) and h14-3-38 (grey) in complex with phosphopeptide (represented as sticks in pink-blue) showing that upon 
phosphorylation g14-3-3 helix I repositions closest to helix H but still doesn't overlap with h14-3-3£ helix I. The h14-3-3£ residues involved in the 
interaction with the peptide are represented as stick (grey). 
doi:1 0.1 371/journal.pone.0092902.g006 



visible suggesting the occurrence of higher oligomers that, 
although, could not be resolved in the gel matrix. In contrast, 
the formation of trimers and tetramers was disfavored in the 
T208A mutant, as well as in polyGlO mutant and, even more, in 
polyG20 mutant (Fig. 7B, left panel), thus suggesting that in these 
mutants C-terminal dimerization could be hampered. Further- 
more, when the incubation and the cross-linking were performed 
in presence of DTT no differences in oligomerization could be 
observed (Fig. 7B, right panel), thus indicating that disulfide 
bridges, if any, are not involved in gl4-3-3 oligomerization. 

Similarly, in native electrophoresis both hi 4-3-3^ and endog- 
enous gl4-3-3 existed exclusively as N-terminal dimers (Fig. 7C). 
Compatibly with the gel resolution limits, for the recombinant 
gl4-3-3, and each of the T214E and R200K mutants, a tetrameric 
form was favored over the dimeric one, with the R200K mutant 
nearly completely tetrameric, as supported by densitometric 
analysis of the protein bands (Fig. 7C-7D). On the opposite, in 
the polyG20, polyGlO and T208A mutants, the dimeric form was 
prevalent, with polyG20>polyG10=T208A (Fig. 7C-7D), further 
supporting that C-terminal dimerization was disfavored. For all 
the tested proteins, binding to target did not affect neither the 



dimerization nor the oligomerization as it resulted by the 
incubation of the proteins alternatively with the Raflp phospho- 
peptide, reproducing a gl4-3-3 consensus mode-1 binding 
sequence of the Raf-1 protein kinase [15], or the corresponding 
non phosphorylated peptide (Rafl) (Fig.7B-7C). However, this 
observation cannot be extended to the endogenous g 14-3-3 that, 
due to the purification procedure, was already in complex with the 
Raflp (Fig. 7B-7C). 

These results strongly indicate that: i) the addition of a C- 
terminal polyglycine stretch, but not the presence of a negative 
charge at position T214, primarily impairs oligomerization of 
recombinant gl4-3-3; ii) mutations R200K and T208A indeed 
affect oligomerization but in opposite directions, i.e. R200K favors 
it whereas T208A does not. 

Finally, to further prove that the observed gl4-3-3 filaments 
were not a crystaUographic artifact and that recombinant gl4-3-3 
can form multimers in solution, the protein was studied with 
transmission electron microscopy. As shown (Fig. 8A), filaments of 
approximately 6 nm of thickness and ranging in length from 0,8 
up to 1,6 |im were visible in the gl 4-3-3 samples. Intersecting 
filaments of greater length (up to 4 |im) were observed when the 
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Figure 7. Assessment of g14-3-3 multimerization in vitro. A) Far-UV CD spectra of 14-3-3 proteins. WT (black), T214E mutant (green), T208A 
mutant (cyan), R200K mutant (violet), polyglycinated 14-3-3 (red). CD spectra were measured with a Jasco J-715 spectropolarimeter (Jasco Ltd, 
Hachioji City, Tokyo, Japan) in 1.0 cm quartz cuvettes between 260 and 195 nm. Samples were 0.1 \ig/^l\ in 10 mlVl Tris/HCI buffer (pH 7.0). B) Cross- 
linking assay. Proteins, pre-incubated with Raflp phosphopeptide, were incubated with DIVIP (upper panels), in presence or absence of DTT, or 
without DIVIP (lower panel), separated on 4-12% reducing SDS-PAGE and immunoblotted with the indicated antibodies. Asterisks indicate: protein 
monomer (*), dimer (**), trimer (***) and tetramer (♦♦**). Molecular size markers (kDa) are on the left. C) Coomassie staining of 12% basic native PAGE. 
Proteins were either pre-incubated with Raf 1 p phosphopeptide or the unphosphorylated peptide Rafi . NativeMark (Invitrogen) size markers (kDa) are 
on the left. Asterisks indicate g14-3-3 and mutants: dimer (**) and tetramer (*♦**). Empty dots indicate h14-3-3C^ dimer ("'). It is to note that, despite 
the theoretical molecular weight of g14-3-3 and h14-3-3iJ are comparable (28.9 and 28.1 kDa, respectively) their theoretical isoelectric point (Ip) are 
much divergent (5.09 and 4.73, respectively) thus contributing to the faster migration of h14-3-3JJ homodimer in native PAGE. Similarly, the presence 
of either phosphorylated T214 in the endogenous g14-3-3 or the glutamic acid at the same position in the T214E mutant decrease the theoretical Ip 
from 5.09 to 5.03 and 5.01, respectively, thus accounting for the increased migration of the proteins in the gel. D) Densitometric analysis of the native 
PAGE presented in (C) performed with ImageJ software. In the graphs: the optic density (OD) of each protein band for each lane is reported as relative 
percentage (%) of the total of optic density (OD) per lane and expressed as relative distribution of g14-3-3 dimer vs tetramer. 
doi:1 0.1 371 /journal.pone.0092902.g007 



protein concentration was increased from 0.1 to 1 |ig/ |il (Fig. 8C, 
8D and 8E). In agreement with native gel and cross-linking data, 
no filament could be observed in the polyG20 sample (Fig. 8B). 
The size of the observed gl4-3-3 filaments was also in good 
agreement with measurement obtained from the solved crystal 
structure (Fig. 8F). Altogether, these data are in favor of the notion 
that oligomerization of gl4-3-3 end up in filamentous structures, 
and non in globular structures, that according to the crystallo- 
graphic data likely depends on domain swapping-mediated C- 
terminal dimerization. 

Discussion 

We report the X-ray structure oiGiardia duodenalis 14-3-3 (gl4-3- 
3). A structural analysis confirmed that the protein has many of the 
features common to the other 14-3-3 family members and also 
highlights functional/structural roles for two peculiar gl4-3-3 



post-translational modifications, i.e. phosphorylation of Thr214 
and polyglycylation of Glu246. 

Despite the quite large number of 14-3-3 crystal structures 
available, we provide the first evidence that a Giardia member of 
the 14-3-3 family, deprived of post-translational modifications, 
forms oligomers in vitro. Analyses of the crystal structure revealed 
that recombinant unmodified gl 4-3-3 forms filaments likely by 
concatenation of the functional dimer. Such process, which was 
further confirmed by transmission electron microscopy, is enabled 
by the combination of N-terminal dimerization and the swapping 
of two hehces between two adjacent g 14-3-3 monomers with the 
formation of a four-helix bundle. Literestingly, this assembly 
retains the peptide-binding groove, but requires the rearrange- 
ment of the HI-loop with respect to other known 14-3-3 structures. 

Domain swapping is dependent on the so called "hinge loop", 
defined as a region of the protein that folds back on itself in the 
monomer but adopts an extended conformation in the domain- 
swapped form. Mutations in this region can in principle modulate 
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Figure 8. Recombinant g14-3-3 form filaments in a concentration-dependent manner. A) Recombinant g14-3-3 (0.1 ^g/|xl) forms 
numerous filaments with moderate length (between 0,8-1 ,6 i^m) and an average diameter of 5,2±0,7 nm, as measured by ITEM soft imaging system 
(OLYMPUS). B) At the same concentration (0.1 |xg/|il), the presence of the polyglycine stretch in the polyG20 mutant prevents filaments formation 
and the protein is found only as amorphous aggregates. C) Longer filaments (more than 3-4 \im) are formed by recombinant g14-3-3 at higher 
concentration (1 |xg/|il). D) and E) are sections of panel C at high magnification (100.000X). Representative fields on the grids are depicted. F) 
Structure-based views of a g14-3-3 filament of concatenated 10 dimers. Dimension (nm) of the filaments are reported. 
doi:1 0.1 371 /journal.pone.0092902.g008 



the propensity of the protein to domain-swap (44). Consistently 
with this notion and our MD simulations, the gl 4-3-3 HI region 
seems to represent the hinge loop, since it hosts residues 
responsible for the "extended" C-terminal conformation, which 
are unique among 14-3-3 family members, such as Thr208 in the 
beginning of loop HI, and Arg200, in helix H. Indeed, mutation of 
Thr208 to an evolutionary conserved alanine and Arg200 to an 
evolutionary conserved lysine affects oligomerization but in 
opposite directions. As expected, T208A slighdy decreases 
oligomerization likely due to the inability of alanine to form 
hydrogen bonds with Gln226 and Arg229 thus weakening the 
interactions between helices H and I. In contrast, the increased 
oligomerization due to the R200K mutation can be explained by 
the ability of the arginine to establish an intra-helix interaction 
with Asn233 and an inter-helix interaction with Asp204 (in helix 
H), thus favoring an equilibrium, whereas the lysine can only 
interact with Asn233 on the same helix (Fig. 3D). 

The exclusively dimeric nature of the in vim polyglycylated g 1 4- 
3-3 and the hampered oligomerization and filament formation in 
both polyGlO and polyG20 mutants strongly suggest that the 
physiological role of polyglycylation is to prevent oligomerization 
of the protein in vivo. This is consistent with the observation that 



gl4-3-3 is always polyglycylated in vivo. In fact, a very short 
polyglycine chain (less than ten residues) can be detectable at the 
gl4-3-3 C-terniinus even in parasite lines over-expressing the two 
gl4-3-3-specific deglycylases [19], supporting a detrimental effect 
of the lack of the polyglycine chain in vivo. This role of 
polyglycylation is in addition to the previously described regulation 
of stage-dependent re-distribution of gl4-3-3 between nuclei and 
cytoplasm already assigned to this PTM [15,18]. The mechanism 
by which polyglycylation impairs oligomerization may be due to 
steric hindrance and/ or disorder promotion. In this regard, the C- 
terminal tail of all known 14-3-3s is enriched in disorder- 
promoting residues (i.e. glycine). In fact, deletion of the last 15 
C-terminal residues of the hi 4-3-3^ induces aggregation of the 
protein in solution at physiological temperature, even if the protein 
retains its ability to bind phosphorylated ligands [45]. Since the 
gl4-3-3 displays a shorter C-terminal tail, polyglycylation may 
introduce disorder in the protein [46-47]. Since there is evidence 
that polyglycine stretches can form extended structures both in 
crystals and in solution, we can speculate that an extended 
polyglycine chain could lay between helices H and I of each gl4-3- 
3 monomer and act as a barrier hindering domain swapping. 
Intriguingly, in the C-terminal part of both the 5. cerevisiae 14-3-3 
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isoforms, BMHl and BMH2, is present a polyglutamine (polyQ) 
stretch (10 and 17 glutamine residues, respectively) and, when 
looking at protein alignment, the polyQ, stretch of both BMHs 
starts approximately in correspondence to the polyglycylation site, 
Glu246, of gl4-3-3 (Fig. 2) [15,18]. As support to our proposed 
role for polyglycylation, it has been shown that the polyQ rich C- 
terminus of BMHs doesn't induce protein aggregation but adopts 
an open extended conformation, is responsible to increase the 
hydrodynamic radius of the BMH molecules, if compared with 
h 14-3-3!^, and it does not lay in the binding groove and cannot 
function as an autoinhibitory binding domain, as described for 
mammalian and plant isoform [48] . In previous works we showed 
that the T2 1 4E mutant displays an increased binding efficiency in 
vitro toward a subset of mode-1 phosphopeptide, if compared to 
the non-phosphorylated recombinant protein [18]. MD simula- 
tions show that the C-terminal tail of the phosphor^'lated protein 
displays a more rigid and "closed" conformation than the non 
phosphorylated one, but yet not sufficient to prevent oligomeri- 
zation. Consistently, the T214E mutant, in which the side chain of 
the glutamic acid mimics the phosphate moiety of the phosphor- 
ylated threonine, is still able to oligomerize. This suggests that the 
primary role of the constitutive gl4-3-3 phosphorylation in vivo is 
rather target binding. In MD simulations, upon phosphorylation 
of Thr214, helix I moves closer to helix H and causes Trp235, 
Asn231 and Leu223 to steadily face the phosphopeptide binding 
groove. Interestingly, these three residues correspond to Trp228, 
Asn231 and Leu216 (hl4-3-3s numbering), which in plant and 
mammalian 14-3-3s are directiy involved in the binding to target 
peptides both in phosphorylation-dependent and -independent 
manner [8]. In contrast, in the non phosphor\'lated gl 4-3-3, the 
interaction of Trp235, Asn231 and Leu223 with the target 
proteins would first require a large conformational change re- 
positioning helix I closer to helix H or would need domain 
swapping of helix I' from another gl4-3-3 dimer. This would likely 
entail a weaker initial binding of gl 4-3-3 to its target. Thus, in vivo 
phosphorylation might promote a faster and/ or much more stable 
interaction with at least a pool of client proteins. This may also 
explain the defect in encystation observed when the non- 
phosphorylatable T214A mutant is overexpressed in Oiardia [18]. 

The determination of the gl4-3-3 structure raises new 
interesting questions. Due to the early evolutionary branching of 
Giardia and the long (■\ ()luti()nars' history of this parasite, why have 
divergent residues hc.v.ii fLxed in the gl 4-3-3 leading to conforma- 
tional rearrangement and C-terminal domain swapping ? Why 
have PTMs been selected over amino acid mutations to overcame 
the risk of oligomerization? Completely deglycylated gl4-3-3 has 
not been observed in vivo so far, however, it cannot be excluded 
that C-terminal domain swapping and limited oligomerization 
may occur in vivo possibly to mediate multiprotein complex 
formation. In this regard, the presence of 14-3-3s in protein 
aggregates have been observed in various human neurodegener- 
ative diseases [49], although the biological significance of 14-3-3 in 
this context is still far to be clarified. Furthermore, arachidonic 
acid, known to favor tau and synuclein polymerization, also 
promotes in vitro oligomerization of hi 4-3-3^ [50]. Intriguingly, 
sequence similarity at the C-terminal region, in particular in the 
HI loop, has also been observed between 14-3-3s [51-52], 
including gl 4-3-3 (data not shown) and a-synuclein, which form 
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Figure S2 Tandem MS analysis of the gl4-3-3 phosphorylated 
peptide (202-219). The MS/MS spectrum of the phosphopeptide 
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